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X  higher-order  parabolic  equation  (PE)  based  on  a  Fade  series  ^J^nd  an  elastic  PE  [-10]  are 
applied  to  wave  propagation  in  the  ocean.  In  contrast  to  the  standard  PE  models  of  underwater 
acoustics,  the  higher-order  PE  provides  accurate  solutions  for  problems  involving  arbitrarily  long 
ranges,  propagation  nearly  normal  to  the  preferred  direction,  and  large  variations  in  sound  speed. , 
The  most  important  applications  of  the  higher-order  PE  are  for  problems  involving  elastic  ocean 
bottoms.  A  new  numerical  approach  based  on  centered  differences  is  applied  to  handle  interface 
conditions.  The  accuracy  of  the  clastic  PE  is  demonstrated  with  benchmark  calculations.  The 
clastic  PE  is  applied  to  a  range-dependent  propagation  problem.  /  ■  , 

1.  Introduction 


The  parabolic  equation  (1]  (PE)  is  very  useful  for  wave  propagation  problems  because  it  can 
be  solved  with  an  efficient  marching  algorithm.  However,  PE  solutions  arc  only  approximate  and 
thus  arc  valid  for  only  certain  types  of  propagation  problems.  The  original  narrow-angle  PE, 
which  is  based  on  a  rational  linear  approximation  r  f  a  square  root  function,  is  fairly  accurate  if 
the  propagation  direction  is  limited  to  within  abo-jt  15  degrees  of  the  preferred  direction.  The 
wide-angle  PE,  (2-1)  which  is  based  on  a  rational  I  near  appro.ximation  of  a  square  root  function,  is 
fairly  accurate  for  propagation  angles  up  to  about  10  degrees.  These  standard  versions  of  the  PE 
are  applicable  to  a  large  category  of  oceanic  wave  propagation  problems.  However,  they  can  not 
handle  problems  involving  very  hard  ocean  bottoms  or  propagation  very  close  to  or  very  far  from 
the  sound  source. 


Several  PE  models  have  been  derived  for  wave  propagation  in  elastic  media.jS- 10)  A  few  of  these 
models  wore  implemented  and  produced  promising  r«^iilts.  (5,9)  However,  no  conclusive  evidence 
has  been  provided  that  the  elastic  PE  is  capable  of  .solving  realistic  problems.  PE  models  that  are 
derived  by  assuming  that  energy  propagates  at  speeds  close  to  a  reference  speed  obviously  can  not 
handle  problems  involving  coupling  between  compressional  and  shear  w.aves.  Furthermore,  clastic 
ocean  bottoms  arc  often  very  hard  and  thus  support  very  wide  angle  propagation. 

To  handle  problems  involving  propagation  at  very  wide  angle.s  or  of  a  superposition  of  difTercnt 
wave  types,  a  higher-order  PE  based  on  a  very  accurate  approximation  of  the  square  root  function 
is  required.  Several  generalizations  of  the  standard  PE  models  have  been  invc.stigated.  {11-1.3) 
In  this  paper,  a  Fade  series  approximation  (13)  is  applied  to  derive  both  .icoustic  and  elastic  PE 
models.  Calculations  are  presented  that  demonstrate  that  these  models  arc  v.alid  for  propagation 
nearly  orthogonal  to  the  preferred  direction  and  for  wave  speeds  quite  different  from  the  reference 
speed.  The  result.s  of  Ilcf.  10  are  implemented  and  extended:  the  elastic  PE  is  gener.alized  to 
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higher  order  and  to  media  in  which  the  Lame  constants  and  density  vary  continuously  and  solved 
numerically  with  Galerkin’s  method;  a  new  numerical  approach  is  used  to  handle  the  interface 
conditions  between  layers. 


2.  Standard  PB  approximations 

We  work  in  cylindrical  coordinates  with  z  being  the  depth  below  the  ocean  surface  and  r 
being  the  horizontal  distance  (range)  from  a  time-harmonic  point  source  of  circular  frequency  u. 
For  now,  the  comple.x  wavenumber  K  =  k  +  it)0  |  k  |  and  density  p  are  assumed  to  depend  only 
on  z,  where  k  =  tu/c,/3  is  the  attenuation  in  decibels  per  wavelength,  c  is  the  sound  speed,  and 
1]  =  (‘lOnlogioc)"*.  We  define  the  reference  sound  speed  cq  and  reference  wavenumber  ko  =  tv/co. 
Cylindrical  spreading  is  handled  by  removing  the  factor  r“l  from  all  field  variables. 

For  kr  »  1,  the  complex  pressure  P  satisfies  the  farfield  equation 


d^P  _  IdpdP 

pdz  dz  ^ 


K^P  =  0, 


(2.1) 


which  factors  exactly  to 


dr 


ik^ 


1  + 


A'*  -  *0  +  £2  -  ilf  ^ 


Jbl 

kl 


(2.2) 


lly  allowing  A'  and  p  to  depend  on  r  in  Eq.  (2.2),  which  we  refer  to  as  PE.x.>  a  leading-order  solution 
is  obtained  for  problems  in  which  range-dependence  is  a  perturbation.  (Mj  The  e.xact  solution  of 
*  PE«  can  be  obtained  in  terms  of  outgoing  coupled  modes.  (15) 

Approximate  solutions  of  PE.x.  arc  obtained  by  approximating  the  square  root  in  Eq.  (2.2)  and 
applying  numerical  methods.  The  narrow-angle  PE  is  obtained  using  the  first  term  of  the  Taylor 
scries 


y/r+n  - 


+ 


t2.3) 


The  plane  wave  factor  exp(ii('or)  i.s  removed  from  P,  and  we  as,sume  that  A’  =  ko  to  obt.ain 
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The  terms  on  the  right  side  of  Eq.  (2.-1)  are  the  refraction  term,  which  accounts  for  variations  in 
sound  speed:  the  don.sity  term,  which  accounts  for  density  variations;  the  loss  term,  which  accounts 
for  sediment  attenuation;  and  the  diffraction  term,  which  accounts  tor  the  vcrtic.al  component  of 
propagation. 
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fiorm.'il  mfiilf  ■solution,  .Solid  curvo-  .ifo  (;i)  (he  narrow  anglo  PF.  ^o^lIion  and  (l>)  the  wnlo  angle 
PF  ■'olutinn. 


riie  narrow-angle  PK  is  arnirate  only  for  a  limiti-d  range  .md  only  if  somid  spced  variations  and 
Iiro|/agation  angles  are  small.  Longer  ranges,  larger  sound  speed  variations,  and  wider  propagation 
angles  ran  be  bandied  by  using  belter  approvimations  for  the  srpiare  ri>ot.  rin-  most  obvious 
approarb  is  to  use  more  terms  in  ibe  Taylor  series,  (llj  Since  tlie  laylor  series  diverges  for  |  J"  1>  1. 
bowever,  many  terms  are  required  for  problems  involving  large  differen.es  in  sound  spi-ed  or  very- 
wide-angle  propagation.  Furlbermore.  this  approarb  is  dilbciilt  to  implement  because  r  is  raised 
to  powers. 
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Table  I:  Comparison  of  Taylor  and  Fade  series. 


X 

4-term 

Taylor 

1-tcrm 

Fade 

2-term 

Fade 

3-term 

Fade 

0.25 

1.11801 

1.11765 

1.11803 

1.11803 

1.11803 

0.50 

1.22412 

1.22222 

1.22472 

1.22474 

1.22474 

0.75 

1.31870 

1.31579 

1.32274 

1.32287 

1.32288 

1.00 

1.39844 

I.IOOOO 

1.41379 

1.41420 

1.41421 

1.25 

1.45639 

1.47619 

1.4990-1 

1  49996 

1.50000 

1.50 

1.48193 

1.54545 

1.57931 

1.58105 

1.58114 

1.75 

1.46078 

1.60870 

1.65523 

1.65812 

1.65831 

2.00 

1.37500 

1.72000 

1.7958-1 

1.80221 

1.80278 

2.50 

0.91943 

1.76923 

1.86124 

1.S6994 

1.87083 

2.75 

0.49545 

1.81481 

1.92376 

1.93510 

1.93649 

3.00 

-0.10156 

1.85714 

1.98361 

1.99817 

2.00000 

%/TTi  - 

II 
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whore  n  is  the  number  of  terms  in  the  Fade  series  and 
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is  based  on  the  1-term  Fadd  series.  The  wide-angle  FE  <lrK*s  not  separate  into  terms  corresponding 
to  simple  phy.sical  proce.ssos. 

We  solve  Eejs.  (2.'1)  and  (2.-'t)  numerically  by  first  discretizing  depth  dependence  with  (lalerkin's 
method  a.s  do,scribed  in  the  Appeiidi.x.  This  approach,  which  is  valid  for  piecewise  continuous 
variatioii.s  in  A'  and  p,  is  ea.sier  to  apply  than  finite  dilference  schemes.  The  resulting  sjslem  is 
then  .solved  with  Crank-iVirolson  integration. 


Figuff  2;  Tran.sinission  los.s  at  j  =  2.'>m  for  a  2'>ih  soiirco  in  a  wavcgiiido  with  poifi'clly  rollprling 
boundaries.  The  da.slied  curvo.H  are  the  normal  mode  solution.  The  solid  rnrvcs  are  t  he  I*lv.  Milnlion 
for  (a)  n  =  1,  (b)  n  =  2,  (r)  n  =  .3,  (d)  n  =  1.  (o)  n  =  .%  and  (f)  »  =  G. 
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Wc  consider  NORDA  Test  Problem  3B,  {16j  which  involves  a  250112  source  at  z  =  99.5m  in 
an  ocean  100m  deep  in  which  c  =  1500m/s.  In  the  sediment,  c  =  1590m/s,  p  =  1.2g/cm^,  and 
=  0.5.  We  use  Greene’s  wide-angle  PE  starter  (•!]  as  an  initial  condition  at  r  =  0.  The  narrow- 
angle  and  wide-angle  PE  solutions  are  compared  with  the  normal  mode  solution  in  Figure  1.  The 
narrow-angle  PE  solution  has  a  large  phase  error.  The  wide-angle  PE  solution  is  very  accurate. 

3.  The  higher-order  PE 


Since  the  Fade  scries  is  valid  outside  the  radius  of  convergence  of  the  Taylor  series,  relatively 
few  terms  arc  needed  for  x  =  1.  Wc  illustrate  this  in  Tabic  I.  The  d-term  Taylor  series  is  better 
than  the  1-term  Pade  scries  for  x  <  I,  but  the  1-tcrm  Padc  scries  is  better  for  i  >  1.  The  2-tcrm 
Pade  scries  and  the  d-term  Taylor  scries  arc  both  correct  to  0(i*)  for  small  x.  Yet  the  2-term 
Pade  scries  is  substantially  better  than  the  d-term  Taylor  scries.  The  3-tcrm  Pade  scries  is  fairly 
accurate  well  beyond  the  radius  of  convergence  of  the  Taylor  series  near  z  =  3. 

The  Pade  series  gives  the  higher-order  PE 


which  wc  refer  to  as  PEn-  Equation  (3.1)  can  be  solved  with  the  method  of  alternating  directions. 
This  approach  involves  n  steps  with  stop  j  requiring  the  solution  of  the  equation 


i£fi.  + 

Pdzdz 

4. 


h’-tS'j 


(3.2) 


Equation  (3.2)  is  solved  with  the  approach  used  to  solve  Eq.  (2.S). 

To  illustrate  the  ability  of  PEr,  to  handle  long-range  and  very- wide-angle  propagation,  wc  con¬ 
sider  a  waveguide  2.50m  thick  with  pre.ssure-rele.aso  top  and  bottom  boundaries  in  which  r  = 
1.500m/s.  A  2511z  point  -source  is  placed  at  z  =  25m.  and  we  take  cq  =  1500m/s.  The  eight 
normal  modes  for  this  problem  propagate  at  approximately  7,  1-1,  21,  29,  37,  -16,  57.  and  7d  degrees 
from  horizontal.  PEn  solutions  (initialized  with  the  norm,al  mode  solution  at  r  =  0)  are  compared 
with  the  norm.al  mode  solution  in  Figure  2.  We  observe  that  the  PEn  solutions  break  down  very 
rapidly  with  r  for  small  n.  However,  the  PEe  solution  is  very  accurate  at  r  =  -Ikm. 

Wc  now  con.sidcr  an  e.xample  that  illuslr.atcs  the  application  of  PEn  for  low-frequency  under¬ 
water  acoustic  propagation  in  deep  water.  In  the  water  column,  we  .'ussiime  the  Munk  profile  [17] 


c(z)  =  c,,  {l  -h  ;r  [2i-^  -f  exp  (-2^5^)  -  l]  }  . 


(3.3) 
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Range  (km) 

Figure  3:  Transmission  loss  at  z  =  200m  fnr  a  lOIIz  source  in  deep  water  over  a  hard  acoustic 
bottom.  The  dashed  curve  is  the  PEi  solution.  The  solid  curve  is  the  PEs  solution. 


Range  (km) 


Figure  -1;  Transmission  loss  at  a  =  2.'>m  for  a  .lOHz  source  in  shallow  water.  The  dashed  curve  is 
the  PEj  solution  for  co  =  1.500m/s.  The  solid  curve  is  the  PEis  solution  for  co  =  .300m/s. 
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Figure  5:  Transmission  loss  at  2:  =  2om  for  a  soft  clastic  bottom.  The  clastic  PE  solution  with 
interface  conditions  (solid  curve)  and  FFP  solution  (dashed  curve)  are  in  nearly  perfect  agreement. 
The  acoustic  PE  solution  (broken  curve)  is  included  to  show  the  importance  of  shear  effects. 

wlijre  It  =  0.0071,  =  1500m/s,  z^i,  =  1000m,  and  II  =  1200m.  The  ocean  depth  is  5000m.  In 

th<  sediment,  c  =  ISSOm/s,  p  =  1.5g/cm^  and  ^  =  0.5.  A  lOllz  point  source  is  placed  at  z  = 
20i)m,  and  we  take  cq  =  1.500m/s.  The  homogeneous  half-space  field  [IS)  is  used  to  initialize  the 
field  at  r  =  400m.  The  Lloyd’s  mirror  beams  produced  by  the  source  propagate  at  appro.ximatcly 
II,  34,  and  70  degrees.  PE|  should  accurately  account  for  the  11  and  34  degree  beams  for  well 
beyond  r  =  20km.  From  the  PEi  and  PEs  solutions  appearing  in  Figure  3,  however,  we  observe 
that  PEi  can  not  handle  the  70  degree  beam  which  is  partially  reflected  from  the  ocean  bottom 
and  makes  a  significant  contribution  to  the  field  for  5km  <  r  <15km. 

To  illustrate  a  possible  application  of  PE„  in  elastic  media,  we  consider  a  problem  for  which  c  = 
I500m/s  in  the  water  and  the  ocean  depth  is  200m.  In  the  sediment,  c  =  1700m/s,  p  =  1.5g/cm^, 
and  p  =  0.5.  A  50Hz  point  source  is  placed  at  z  =  2.5m.  The  PE2  solution  for  cq  =  1500m/s  and 
the  PEi3  solution  for  co  =  300  m/s  appear  in  Figure  4.  The  agreement  o'.'  the  solutions  suggests 
that  a  higher  order  elastic  PE  based  on  the  Padd  series  would  handle  both  (ompressional  and  shear 
waves  simultaneou.sly. 


4.  The  elastic  PE 


The  farfield  eciuations  for  the  horizontal  and  vertical  displacement  u  and  r  in  an  ela-slic  medium 
are  [19j 
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where  the  dil'  .tion  A  is  defined  by 


dp  dll  dpdv  _ 

'dz'Oz  dr  ~ 

('1.1) 

dpdv  ^ 
dzdz  ^  dz 

('1.2) 

A  =  +  §.  ('»-3) 

Loss  is  handled  by  using  the  complex  compressional  and  shear  speeds  Cp  =  Re(Cp)/(l  +  iq0p) 
and  c,  =  Re(c,)/(1  +  iq0s),  where  /Jp  and  0,  are  the  compressional  and  shear  attenuation  and 
q  =  (-lOx  logioc)"*.  The  Lame  constants  A  and  p  arc  related  to  the  wave  speeds  and  the  density  p 
by  A  =  p  (c^  -  2c;)  and  ;i  =  pcj. 

We  dilTercntiatc  Eq.  ('Ll)  with  respect  to  r  and  Eq.  ('1.2)  with  respect  to  z  and  sum  them 
using  Eq.  (■4.3)  to  obtain 

(A  +  2/i)^  +  (A  +  2/i)^  +  pu;*A  +  2|^|p  +  + 

Equations  (>*.2)  and  (  l.'l)  provide  a  coupled  system  of  the  form 


which  reduces  to  tlie  formulation  used  in  Ref.  10  for  a  homogeneous  medium.  Since  the  operators 
L  and  M  commute  with  d/Or.  we  may  factor  Eq.  (  l-l)  to  obtain  the  outgoing  wave  equation 

The  plane  wave  factor  exp(jA(,r)  is  removed  from  A  and  t>  and  the  square  root  in  Eq.  (  LG)  is 
approximated  with  the  Fade  series  to  obtain  the  higher-order  elastic  PL 


1  XT  «,.n/'-'('U  -  kjl.) 

~  °  kl  +  b,,„lrHM-kih)  I  I 
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Figure  6:  Transmission  loss  at  z  =  25m  for  a  soft  clastic  bottom.  The  elastic  PE  s  >lution  with 
Galerkin’s  method  for  the  interface  (solid  curve)  and  the  FFP  solution  (dashed  curve)  arc  in  fairly 
good  agreement. 

We  solve  the  clastic  PE  numerically  with  the  method  of  alternating  directions  and  Crank- 
Nicolson  integration  to  march  the  solution  in  r.  The  depth  operators  arc  replaced  with  matrices 
using  Galcrkin's  method  within  layers  in  which  A,/i,  and  p  are  continuous  as  described  in  the 
Appendix.  Conditions  for  the  continuity  of  the  displacements  and  the  normal  and  tangential 
stresses  arc  used  at  interfaces  between  solid  layers.  At  a  fluid-solid  interface,  however,  horizontal 
displacement  is  not  required  to  be  continuous. 

Let  us  consider  a  fluid-solid  interface  at  i  =  z,  (the  grid  points  z,  are  defined  in  the  Appendix) 
and  introduce  the  subscripts  te  and  5  for  the  fluid  (water)  and  solid  (bottom)  layers.  We  use  the 
physical  values  Au(jj),Ai(Zj),  V4(z,).  A,..(Zj_i),A)i(j,+i),  and  «'a(z^+i)  and  the  nonphysical  v.ilues 
^k-(^j+i)»^’‘(*j-i)t  <‘nd  ns(zj.])  in  the  discreti/eil  equations  of  motion  in  both  layers  at  z  =  z,. 
Wo  also  use  the.se  values  to  form  difTerenres  at  z  =  z,  for  the  following  conditions  for  continuity  of 
vertical  displacement,  normal  stress,  and  tangential  stre.ss  (10): 


—  (A„A«.)  -p  =  0 

(■I.S) 

A^  A,^.  =  As  As  +  2n>.~ 

Oz 

(1.9) 

=  0. 

Oz  Oz  Oz 

(1.10) 
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Figure  7:  Transmii.$ion  loss  at  ;  =  25ni  for  a  hard  clastic  bottom.  The  elastic  F'E  (solid  curves) 
and  the  FFP  solution  (dashed  curves),  (a)  The  solutions  do  not  agree  for  n  =  1.  (b)  The  solutions 
are  in  good  agreement  for  n  =  3. 


T'he  nonphysical  values  are  eliminated  in  the  equations  of  motion  at  -  =  using  the  diifercnce 
formu!a.s  for  Eqs.  (-1.8),  (d.O),  and  (<1.10),  A  solid-solid  interface  is  handled  similarly  using  the 
interface  conditions  described  in  Ref.  10. 


In  a  fluid  medium,  P  =  \A  can  be  initialircd  at  r  =  tq  with  the  homogeneous  half-space  field 
[l.'i]  Pi,  for  a  point  source  at  r  =  lo 


Pk  = 


luld. 


1 

—  e.\p 
d+ 


iwd^. 

^ep(*o) 


1 

-exp 


(■1.11) 
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4  =  rg  +  (4.12) 

where  fcro  >>  1.  The  half-space  field  accounts  for  the  direct  arrival  and  the  surface-reflected  arrival. 
This  PE  starter  is  generalized  to  clastic  media  by  applying  the  half-space  quantities 


~  A(zo) 

(4.13) 

1  dPh 

Vh  - - 5  -r- 

dz 

(4.14) 

to  initialize  A  and  v.  This  starter  is  valid  only  if  the  source  is  at  least  a  few  wavelengths  away 
from  the  ocean  bottom.  We  define  transmission  loss  to  be  — 20!og,o  |  A(2o)A  J  -flOlogior.  This 
definition  corresponds  to  the  usual  definition  of  transmission  loss  in  a  homogeneous  fluid. 

We  first  consider  a  problem  involving  a  50Hz  source  at  z  =  25m  in  an  ocean  of  depth  100m  in 
which  Cp  =  1500m/s.  In  the  uniform  clastic  sediment,  Cp  =  1700m/s,  c,  =  800m/s,  p  =  Ig/cm®, 
fip  —  0.5,  and  0,  =  0.5.  For  the  numerical  solution  of  the  clastic  PE,  we  used  the  values  n  =  1, 
Ar  =  5m,  Az  =  0.5,  and  tq  =  100m.  To  eliminate  an  instability  in  the  numerical  solution,  it  was 
necessary  to  extend  the  computational  domain  2km  below  the  ocean  surface  at  which  we  assumed 
that  botli  A  and  v  vanish.  Both  the  comprcssional  and  the  shear  attenuation  were  increased  linearly 
to  10  in  the  lower  100m  of  the  domain  to  prevent  reflections.  Transmission  loss  generated  with  the 
clastic  PE  and  with  an  FFP  model  [20,21]  appears  in  Figure  5.  The  solutions  arc  in  nearly  perfect 
agreement. 

We  also  performed  an  clastic  PE  calculation  for  this  problem  using  the  diflercncc  equations 
obtained  with  Galerkin’s  method  to  model  the  ocean  bottom  interface.  This  approach,  which 
produces  accurate  results  for  acoustics  problems,  [22]  is  easier  to  implement.  Although  L~'  does 
not  exist  where  /i  vanishes,  the  r  solution  of  Eq.  (4.2)  is  also  a  e  solution  of  Eq.  (4.7)  in  a  fluid. 
Thus  the  elastic  PFl  should  be  v’alid  for  p  =  0.  We  observe  from  the  transmission  loss  data  appearing 
in  Figure  G  that  Galerkin’s  method  handles  discontinuities  fairly  well.  This  is  consistent  with  the 
discussion  in  Ref,  23.  We  deduce  that  Galerkin’s  method  should  handle  continuous  variations  very 
well. 

We  now  consider  a  problem  involving  a  20Hz  source  at  z  =  25m  in  an  ocean  of  depth  600m  in 
which  =  I500m/.s.  In  the  uniform  ola.stic  sediment,  r,,  =  3400m/s,  c,  =  1700m/s,  p  =  2g/cm^,idp 
=  0.5,  and  3,  =  0.5,  1  he  hard  ocean  bottom  supports  very  wide  propagation  angles  in  the  water 
column.  For  the  numerical  solution  of  the  elastic  PE,  wo  extended  the  computational  domain  to  z 
=  10km  and  used  the  values  Ar  =  10m,  Az  =  Im,  and  =  200m.  The  attenuation  was  increrised 
linearly  to  100  in  the  lower  1km  of  the  domain.  The  elastic  PE  solution  appears  in  Figure  7  for 
n  =  1  and  for  n  =  3.  The  elastic  PR  can  not  handle  the  wide  propagation  angles  for  n  =  1. 
However,  the  n  =  3  solution  is  in  good  agreement  with  the  FFP  solution. 

In  PE  modeling,  range-dependent  problems  are  handled  approximately  by  allowing  the  en¬ 
tries  of  the  matrices  for  the  depth  operators  to  depend  on  range.  This  approach,  which  is  easily 
implemented,  produces  accurate  results  for  acoustics  problems  if  range  dependence  is  sufficiently 
gradual. 
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Loss  (clB  re  Im)  Loss  (dB  re  Im) 


180 


M.D.  Collins 


Figure  9:  Transmission  loss  at  z  =  30m  for  a  sloping  elastic  bottom.  The  clastic  PE  solution  for  n 
=  10  (solid  curve)  and  n  =;  1  (dashed  curve)  .and  the  acoustic  PE  solution  (broken  curve). 

This  approach  should  also  be  valid  for  weakly  range-dependent  problems  in  clastic  media. 

We  modify  a  range-dependent  benchmark  problem  12'1]  to  illustrate  the  higher-order  clastic  PE 
model  for  a  range-dependent  problem.  A  25112  source  is  located  at  r  =  100m  in  an  ocean  in  which 
c  =  150Cm/s  and  the  depth  d  is  given  by 

In  the  sediment,  Cp  =  1700ni/s,  p  =  1..5g/cm*,  .and  Py  =  0.5  for  the  origin.al  acoustic  problem.  We 
take  c,  =  800m/s  and  >3,  =  0.5.  Due  to  the  sloping  ocean  bottom,  a  significant  amount  of  energy 
is  transmitted  into  both  shear  and  compressional  w.avcs  in  the  bottom.  Thus  one  might  e.xpect  the 
71  =  1  cla,stic  PE  to  have  difficulty  handling  this  problem. 

We  extend  the  grid  down  to  z  =  10km  ,and  take  Ar  =  5m,  Az  =  Im.  Tq  =  lOOm,  and  Cq 
1.500m/N.  Contour  plots  of  transmission  loss  appear  in  Figitre  S  for  77  =  1  and  77  =  10.  The  77  =  1 
solution  breaks  down  in  the  sedime7it  because  it  does  not  handle  the  transmitted  shear  waves.  The 
77  =  10  solution  does  not  appear  to  break  <lowii  in  the  sediment.  Transnrission  loss  curves  appear  iit 
Figure  9.  We  observe  that  the  ehuslic  PE  solution  exhibits  more  loss  than  the  acoustic  PE  solution 
and  that  the  n  =  1  .and  n  =  10  solutions  are  in  good  ;igri‘ement  in  the  water  rolurnn. 
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5.  Conclusion 


Higher-order  acoustic  and  clastic  PE  models  based  on  a  Fade  scries  have  been  derived  and 
solved  numerically.  The  higher-order  PE  models  produce  accurate  solutions  for  propagation  nearly 
orthogonal  to  the  preferred  direction  and  for  large  variations  in  sound  speed.  In  particular,  the 
sound  speed  may  deviate  far  from  the  reference  sound  speed.  The  clastic  PE  allows  hard  ocean 
bottoms  as  well  as  continuous  variations  in  the  Lame  constants  and  density.  It  appears  to  handle 
propagation  of  a  superposition  of  different  types  of  clastic  waves.  Galerkin’s  method  was  used 
to  discretize  the  clastic  PE,  which  has  relatively  complicated  depth  operators.  A  new  numerical 
approach  was  used  for  the  interface  conditions  between  layers.  The  half-space  PE  starter  has  been 
generalized  to  clastic  media.  Benchmark  calculations  demonstrate  that  the  clastic  PE  and  the  half- 
space  elastic  PE  starter  produce  accurate  solutions  for  range-independent  problems.  The  elastic 
PE  was  applied  to  a  range-dependent  problem. 

The  half-space  starter  is  not  valid  for  very  low  frequencies  or  for  sources  close  to  the  ocean 
bottom.  Since  these  kinds  of  problems  involve  strong  coupling  into  shear  and  interface  waves,  an 
clastic  PE  starter  that  handles  them  would  be  useful.  With  the  boundary  conditions  we  used  at 
the  bottom  of  the  grid,  the  clastic  PE  is  stable  only  if  the  domain  is  truncated  at  a  largo  depth. 
Bottom  boundary  conditions  that  allow  smaller  truncation  depths  would  improve  the  efficiency  of 
the  model.  It  might  be  possible  to  accurately  handle  interfaces  with  Galcrkin’s  method  by  slightly 
modifying  the  formulation  of  the  elastic  PE  and/or  by  using  dilTerent  basis  functions.  This  would 
improve  the  simplicity  of  the  model. 
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